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Abstract 

The relation between collective modes and the phase transition in 
low density nuclear matter is examined. The dispersion relations for 
collective modes in a linear approach are evaluated within a Landau- 
Fermi liquid scheme by assuming coexisting phases in thermodynam- 
ical equilibrium. Temperature and isospin composition are taken as 
relevant parameters. The in-medium nuclear interaction is taken from 
a recently proposed density functional model. We found significative 
modifications in the energy spectrum, within certain range of temper- 
atures and isospin asymmetry, due to the separation of matter into 
independent phases. We conclude that detailed calculations should 
not neglect this effect. 
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1 Introduction 



A known feature of the low density regime of the nuclear interaction is that 
matter undergoes a phase transition of the liquid- gas type. There is significa- 
tive evidence of its existence in nuclear collisions experiments at medium and 
high energies. It would also manifest in the crust of neutrons stars, where 
nuclear correlations have a determinant effect on the transport properties 
and, therefore, in the cooling process of stellar matter. 
Under the conditions of interest, namely low density and temperature, matter 
is essentially composed of protons, neutrons, and circumstantially leptons. 
Therefore, we have a strongly interacting binary system, with more than one 
conserved charge. A characteristic aspect of this kind of phase transitions is 
that conserved charges do not distribute homogeneously among the coexist- 
ing phases. This property give rise to the isospin fractionation observed in 
heavy ion collisions. 

Another interesting phenomenon of the bulk nuclear matter is the exci- 
tation of collective modes with a definite energy spectrum. This is a well 
studied subject because of its multiple consequences. For example, in astro- 
physical applications it has been established that collective modes propagat- 
ing in neutron matter modify substantially the neutrino scattering rates pQ. 
A close relation exists between collective modes and the liquid-gas nuclear 
transition, in particular the unstable ones are the precursor of the thermo- 
dynamical instability. There is a great amount of research concerning on one 
hand the collective modes and on the other one the low density phase tran- 
sition in nuclear matter, however a few number of investigations deal with 
the relation between them [2j|3]. For instance, in [3] the spinodal instability 
and the collective modes are studied in a neutral mixture of nucleons and 
electrons, including both iso-scalar and iso-vector density fluctuations. This 
is an appropriate scenario if the system evolves out of equilibrium. However 
for characteristic times large enough, a succession of coexisting phases must 
be taken into account. 

In this letter we aim to study the propagation of collective modes in a medium 
composed of different phases in thermodynamical equilibrium. In contrast to 
[21 13] we have chosen the coexistent phases as the reference state over which 
perturbations are applied. We assume the system is well described in a Fermi 
liquid scheme, and a linear approximation to the transport equation is taken. 
For this purpose we used a model of the nuclear interaction recently proposed 
[I], inspired by the Kohn-Sham density functional approach. Its concise for- 
mulation enables us to include temperature and isospin composition, and to 
concentrate on the physical aspects instead of calculational complexities. 
Our conclusions could be related to the evaluation of cooling rates mediated 
by neutrinos in proto-neutron stars. 



2 The model and its equation of state 



The Density Functional Theory (DFT) seems to be an appropriate field 
where different formulations of the in-medium nuclear interaction could con- 
verge. Dissimilar approaches like Brueckner-Hartree-Fock calculations with 
free-space two nucleon potentials, relativistic field models in mean field ap- 
proximation, non-relativistic effective forces or in-medium chiral perturbation 
theory, provide material for DFT. It has the advantage to yield accurate re- 
sults within relatively simple calculations. 

The Kohn-Sham DFT reduces the description of a complex interacting sys- 
tem to a simple energy functional resembling that of independent particles 
moving in an external potential. This procedure is justified by the Hohen- 
berg and Kohn theorem [5]. 

Within this framework, it has been proposed recently an energy functional 
for finite nuclei jl], which can be decomposed according to E = K + Eso + 
E!j£ t + E-£j! + E c . The first and second terms contain the uncorrelated con- 
tributions to the kinetic energy and the spin-orbit splitting respectively, E c 
stands for the Coulomb energy, and E°^ t + E^f represents the independent 
nucleons moving in a mean nuclear potential. While E°£ t describes the bulk 
matter behavior, the remaining term collects finite range effects. Since we 
are concerned with infinite homogeneous matter, only kinetic and bulk terms 
are retained, in particular we have 



ETm = J d 3 r [P.(x)(l - w 2 ) + P n (x)w 2 ] 



it 



Here n = rii + n 2 is the total baryonic density number, sum of the proton rii 
and neutron n 2 densities, w = (n 2 — ni)/n is the isospin asymmetry fraction, 
and P s {x), P n {x) are interpolating polynomials for symmetric and neutron 
matter respectively. They are written in terms of x = n/riQ with no = 0.17 
fm~ 3 , the number density at saturation, 

p ( r ) = / ^ ) x - 1 m 

aK J I P s (l) + a!(a:-l) + a 2 (a;-l) 2 , x>\ K> 

Pn(x) = j^b^X k (2) 
1 

The coefficients ai, and a 2 can be consulted in [4]. The validity 

of this expansion extends up to x ~ 1.4. 

The polynomials have been obtained by adjusting the correlation part of the 
energy per particle for symmetric and pure neutron matter, and then making 
a quadratic approximation in w. A detailed description of the formalism can 
be found in [6]. 

Within a Fermi liquid approach, the single particle spectra can be found by a 
functional derivative e&(p) = SE/Sf®, with f®(T,p) = [1 + exp(e;,(p) — /i^/Tp 1 
the equilibrium statistical distribution function for a state of isospin 6 = 1 



(protons) or b = 2 (neutrons) 



e&(p) 



2m 



+ [P s (x) + xP' s (x)\ (1 - w 2 ) + [P n (x) + xP' n {x)\ w 2 + 
2w(w + I b ) [P s (x)-P n (x)}, 



where m is the degenerate nucleon mass, and If, = 1(— 1) for protons (neu- 
trons). The prime symbol indicates derivation with respect to x. 

For a given temperature and particle density the corresponding chem- 
ical potential is found by the relation: 



The thermodynamical pressure is given by: P = Y^ib^bi^b — £ + TS, where 
£ — E/V is the energy density, and the entropy per unit volume is 



The numerical factor 2 takes into account the spin degeneracy. 
A homogeneous system at temperature T and isospin composition ni, n 2 
remains thermodynamically stable if the free energy per unit volume J 7 is 
lower than any linear combination of energies corresponding to independent 
phases, satisfying the conservation laws, i. e. 



F{T, n u n 2 ) < A JF{T, n' x , n 2 ) + (1 - A) JF{T, n'[, n 2 ), < A < 1, (3) 



where rik = A n' k + (1 — A) n'l, (k = 1,2), set up the conservation requirement 
for particle and isospin number [7]. 

Alternatively, the condition ([3]) can be stated as a set of two equations [8], 
DetJ > 0, Tr# > 0, with fo? = d 2 J- '/dntdrij. They determine the spinodal 
region in the phase diagram. 

In a phase transition the emerging phases satisfy the Gibbs condition for equi- 
librium coexistence: n' x {T, n'un^) = fi'{(T,n",n,2), fJ, r 2 (T, n^, n' 2 ) = ^{T, n'[, n'^), 
PiTjn'^n^ = P{T, n", n'^), which fix the boundary of the binodal region. 
The binodal encloses the spinodal region since phase separation can occur 
while the condition ((3]) is locally satisfied |7]. For a given temperature, the 
coexistence conditions can be easily interpreted in terms of isobaric repre- 
sentations of the proton and neutron chemical potentials as functions of the 
proton abundance Y = ni/n, as shown in Fig. 1. The points lying on these 
curves which are the vertices of a rectangle correspond to those equilibrium 
states which coexist in a phase separation. From these four points, the pair 
of neutron-proton chemical potentials on the left (right) correspond to the 
low (high) density coexisting phase. As it is explained below, these points 
are used to determine the equation of state within the binodal zone. 
We have studied the equation of state of the system with variable isospin 
fraction at low temperatures, some results are shown in Fig. 2. We have 





used the Gibbs prescription and the conservation laws to evaluate the two- 
phase coexisting region, which produces an almost linear behavior between 
the end points corresponding to a single phase state. The free energy of the 
system is given by J-(T, n 1; n 2 ) = A J-"(T, n[, n' 2 ) + (1 — A) J-(T, n", n'^), where 
A, n[, n' 2 , n", n' 2 ' had been determined previously by the Gibbs construction. 
It must be pointed out that for symmetric matter each phase keeps w = 
with a constant partial density and, as the total density increases, only the 
relative volume fraction changes. That is, the binodal construction reduces 
to the standard Maxwell procedure. In the upper panel of this figure the 
pressure as a function of the baryonic density is shown for several isospin 
asymmetries and T=2 MeV. The system becomes unstable at very low den- 
sity, showing a negative compressibility. The only exception corresponds to 
pure neutron matter. In each case the dashed line represents the theoretical 
prediction without phase separation. In the lower panel the free energy for 
T=2 MeV is shown for a selected set of asymmetry values. As in the pre- 
vious case, dashed lines correspond to the unphysical situation. It can be 
appreciated that the phase separation effectively causes a lowering of T . 

For a given temperature T the pair of states satisfying the Gibbs condition 
determine a closed curve as the pressure ranges within the set < P < P C {T). 
This curve is the intersection of an isothermal with the binodal surface, which 
extends from zero to the critical temperature T c . We exhibit in Fig. 3 projec- 
tions of the binodal into the P — Y plane, corresponding to T = 2, 6, and 10 
MeV. The area enclosed by the curve decreases with temperature, reducing 
to a point at T c . In our calculations we have found T c ~ 12.8 MeV. The 
full curve can be separated into two sections with common end points, one 
having Y = 0.5 and the other one with Y = Y c , the proton abundance at the 
critical pressure P C (T). Most of the coexisting phases at a given temperature 
are represented by a pair of points, located in different sections. In such a 
case lower Y values are associated with the less dense phase. This behavior is 
compatible with the isospin distillation observed in multifragmentation: the 
most dense phase approaches to isospin symmetric matter while the other 
one keeps the neutron excess. 

For further development, we describe here the Landau parameters of the 
model. They are defined in terms of a second variation of the energy density 
[9], namely F BB i = V 2 d 2 £ /dfs dfs> ■ In order to simplify the notation we 
have resumed in only one symbol the indices for isospin, spin and linear 
momentum, i. e. B = (b, s,p). The Landau's parameters are defined as the 
Fourier coefficients of an expansion in terms of the Legendre polynomials 



where v — p • p'/{pp')- 

Within the energy functional model used here, £ does not depend on z/, 




therefore the non-zero components correspond only to I = 0, giving 

2 

^BB' - ^66' = _ [-^0 + (— 1) Gbv\: 

F Q = w 2 (P n - P s ) + {l + w 2 )xP' s -w 2 xP^ 

G bV = {P n -P s ) + 2 5 w w[P n -P s + x{P' n -P' s )} (4) 



3 Collective modes at finite temperature 

Collective modes are associated to local density fluctuations that propagate 
in the nuclear mean field. These fluctuations are the effect of small perturba- 
tions of the occupation distribution f B around its equilibrium form f B . As a 
conserved charge the particle density n b can be written as a summation over 
the level occupation using either f B or f B , 

^ = ^E^ T ) = ^E^cn (5) 

s,p s, p 

The momentum can be regarded as discrete by imposing appropriate bound- 
ary conditions. 

Local density fluctuations cause small deviations from the equilibrium 
distribution, this perturbation is assumed of periodic oscillatory nature 

f B (T,v,t) = f B (T) + Sf B (T,v,t) 

= /o (T ) + fsiT) [I- f B (T)} UBe i (tl .r-„t) (6) 

Our present ansatz for the perturbation 6f B (T, r, t) ensures that in the limit 
T — > only the levels around the Fermi surface contribute to the zero mode 
propagation. The quasiparticle energies change accordingly, since in the 
Landau-Fermi liquid model they are obtained in a self-consistent manner 
from the distribution function j£l [10] 

6e B (T, r, t) = \- ^ F B B , 5f B > (T, r, t) (7) 



B' 



The propagation of these perturbations at low temperature is governed by 
the Landau's kinetic equation. For temperatures well below the Fermi energy 
of the system T <^ €f, the collision term can be neglected [IT] 



O/b dfs_ des_ _ dfs_ de^ = Q 
dt dr dp dp dr 

Introducing Eq.Q into Eq.Q and keeping only linear terms in the fluctua- 
tions, we obtain 

dSfs 86 f B de B _ &f% 85 e B = Q 
dt dr dp dp dr 

which can be further reduced to 

(p-q)v-r. fS-m [i - f%.(T)} 



m 



u B - — y F B b> ~ u B > = 0. (10) 



Since the Fb b 1 do not depend on p, see Eqs. (jlj), we can write 
where 



B> 



T 



u B > 



W b > 



w b 



ir 2 T 



dpp 2 f b °(T,p) [l-f b °(T,p))u b (p) 



(11) 



(12) 



and we have used ub = u b (p), = f b (T,p) to show explicitly their 

dependence on p. Therefore the system fflOl) can be rewritten in the form 



w b + C b {T) F bv wv =0 (6=1,2) 



(13) 



b'=i 



with 



cm 



i 



vr 2 T 



fioo(s) 



dpp 2 f»(T,p) [l-f»(T lP )}n 00 (V z /v p ) (14) 
where v p = p/m, V z = oj/q is the phase speed of the collective mode, and 

's-r 

= I In 1 

i-i 2 y-s 



s + l 



(15) 



is the Lindhard function [9]. 

We shall consider the possibility of damped waves, in which case the real zero- 
mode frequency u acquires an imaginary component, namely u — > oj — %r\. 
Introducing 8 = rj/u, in the case of slightly damped modes where < 5< 1, 
we have to leading order 



ReC b 
lmC b 



* 2 T J 

5 



+ 



7T 2 T 

sgn(J) 
2tt 



dpp 2 f b °(T,p) [l-A°(T,p)] 
d P p 2 f b °(T,p) [l-f b °(T,p)} 
m 2 V z f b °(T,mV z ) 



1 + 

2v„ 



v z 


- In 


v z - 


- Vp 


2% 


V z + v p 


In 


V z - 


- v p 




V z + v p 


+ 



V 2 



V 2 -v 2 
p 



(16) 



On the other hand, instability modes are characterized by u = and rj < 0, 
in which case 



Re C b 



ir 2 T 



dpp 2 f b °(T,p) [l-f b °(T,p)] 



1 - 



rj_ 

qv p 



arctan (qv p /rj) 



Im Ct = 



(17) 



The proper frequencies are identified with the roots of the determinant of the 
system of equations ({TBI , which, neglecting higher orders in Im C b , reduces 
to 



[1 + F^ Re d] [1 + F 2 U 2 Re C 2 \ - F^ 2 Re C x Re C 







F X1 [1 + F 2 ° 2 Re C 2 ] Im d + F 2 ° 2 [1 + F 1X Re &] Im C 2 

-F° a 2 Im {Ct C 2 ) = (18) 



It is easy to show from Eqs. ljTTj) and (fl2j) . that in the present model the 
rate of proton to neutron amplitudes is momentum independent, namely 
U1/U2 = W1C2/W2C1. As done at zero temperature [10], the isospin char- 
acter of the proper modes can be classified as iso-scalar (iso-vector) for 
Re(u 1 /u 2 ) > (< 0), respectively [T2] , 

We have verified that along the phase transition the system is composed of 
two stable independent phases. In particular the instability modes have com- 
pletely disappeared and only stable eigenmodes are present. Furthermore, 
in all the cases considered here, no zero sound modes are found propagating 
in the lower density phase of the coexisting region. The reason is that -F 2 ° 2 
takes always negative values in this phase. 
It is convenient to define an average Fermi velocity Vp 

v 1 ^n b [f~dpp*f b (T,p)[l-f b (T,p)] 
F m^n\j f °°dpp*f°(T,p)[l-f°(T,p)) 1 ] 

as a reference to estimate deviations from the Fermi surface, which gives a 
measure of the validity of the approximations. 

In Fig. 4 we display the typical low temperature dispersion relation for dif- 
ferent isospin composition. At finite temperature the proper collective modes 
arise in pairs. One branch is slightly damped, meanwhile the other one prop- 
agates without dissipation. 

For the sake of comparison, we also include the unphysical results obtained 
by neglecting the phase separation. This situation give rise to, among others, 
the instability modes. In such a case we plot \rj\/q instead of V z . 
For isospin symmetric matter (Fig. 4, w = 0) the mixed phase extends up 
to near n , and as it was noticed each phase stay at constant partial density 
during the transition. As a consequence the collective modes propagate at 
constant speed V z . At low densities there are two branches of iso-vector char- 
acter, which propagate in the liquid phase and continue for higher densities, 
beyond the transition. In addition, a bivaluated stable iso-scalar mode ap- 
pears at densities about Lino- For iso-vector as well as for iso-scalar modes, 
the branch above Vp corresponds to undamped motion. The dissipation in 
the iso-vector wave has an average value of 5 ~ 0.15, meanwhile the iso- 
scalar one has 5 ~ 0.09 — 0.19. Therefore, damping remains very small in 
agreement with our assumptions. 

If phase separation in symmetric matter is disregarded (dashed lines), a dou- 
ble iso-vector mode appears at very low density with small values of V z . It 
grows monotonously with density and joins smoothly with the solutions for 
the mixed stable phase. It has a moderate damping 5 ~ 1 at low densi- 
ties, decreasing to 5 ~ 0.15 near n^. At the same time, a pair of iso-scalar 
excitations arise at densities below 0.8 no- The lower one corresponds to a 
instability mode reflecting the existence of the spinodal region. In general 
the instability modes evolving out of thermodynamical equilibrium are iso- 
scalar since both isospin components are equally affected, in a process which 
eventually leads to clusterization [12J. The damped stable scalar mode has 
a dissipation coefficient growing from 5 ~ 1 to 4 as the density is increased, 



mainly because it approaches the unstable configuration. 
Turning to asymmetric matter (Fig. 4, w — 0.2 — 0.6), this scene changes 
gradually as the density range of coexistence reduces with increasing w. At 
low and medium densities only iso-vector modes propagate with almost con- 
stant velocity V z ~ 0.3. For w > 0.5 these excitations are strongly suppressed 
around the medium density zone. In fact in the denser coexisting phase, the 
only one which can sustain collective motion, the neutron-neutron Landau 
parameter F® 2 decreases and can take negative values when w increases. 
For example, in the case of w = 0.6 the parameter F 22 vanishes around 
n pa 0.53no, and remains negative up to n ~ 0.81no. This explains the sup- 
pression of the density collective modes at medium densities. For supranor- 
mal densities the parameter F 22 grows again, and four collective zero modes 
reappear. For high isospin asymmetry these four branches are of mixed 
character, that is, they change from iso-vector to iso-scalar as the density 
increases. 

It must be pointed out that in neutron rich environments such as w = 0.8 
(not shown here), two iso-vector modes with V z ~ 0.3 persist in the very 
low density regime. This scenario of collective waves propagating in highly 
asymmetric nuclear matter is expected to hold within the inner crust of neu- 
tron stars. The precise structure of this crust is still doubtful, but condensed 
nuclear droplets immersed in a uniform fluid environment of almost pure neu- 
tron matter constitutes a plausible assumption [T3J[T3]. Although Coulomb 
effects are responsible of the existence of the droplets, the mean field prop- 
erties of the asymmetrical extended fluid are governed by the nuclear forces. 
In fact, a liquid-gas coexistence of this low density asymmetric nuclear mat- 
ter environment prevents further instabilities, as it was previously stressed. 
Moreover, this fluid phase can sustain coherent density fluctuations, which 
open a channel for neutrino dispersion. These collective excitations are simi- 
lar to the spin- wave propagation investigated in pQ. Therefore it is reasonable 
to expect a strengthening of the effects discussed there, in particular a more 
pronounced reduction of the in-medium mean free path of neutrinos. A fur- 
ther investigation will be presented elsewhere |15j . 

If the phase separation is not taken into account (dashed lines), a stronger 
suppression of the collective modes is found, mainly because of the more 
pronounced decrease of F 22 with asymmetry. In this situation medium and 
high density excitations are affected, and already for w = 0.4 all of them 
have practically been extinguished. Only the unstable and stable iso-scalar 
modes of sub-saturation densities survive in this case. 

In general for each damped mode the parameter 5 increases a bit with grow- 
ing asymmetry w, and this enhancement is more important when only the 
unstable one phase is considered. We can conclude that the binodal phase 
transition tends to stabilize the density collective zero modes at low and 
medium densities. 



4 Conclusions 



We have applied a density functional model of the nuclear interaction fitted to 
describe asymmetric nuclear matter properties, to determine the coexistence 
regime of the liquid-gas phase transition at densities below the saturation 
value no- This study has been performed at several temperatures and isospin 
asymmetries. The binodal region is diminished as the temperature increases, 
till its critical value T c ~ 12.8 MeV. For a given temperature this region also 
decreases with growing isospin asymmetry. 

We have also studied the propagation of zero sound modes at finite tem- 
perature, using a linearized collisionless Landau kinetic equation. As the 
method requires a stable reference state which supports density fluctuations, 
we have chosen the coexisting phases in the binodal instead of the unstable 
spinodal region. We have found that collective modes are supported only 
by the denser phase of the coexistence region, which favors their propaga- 
tion at low total densities. Because of this, the phase speed of the collective 
excitations remains almost constant in a wide density range. Furthermore, 
since in the denser phase the proton to neutron fraction is higher than in the 
lighter one, it favors the stability of the density zero modes in matter with an 
overall neutron excess. This could have significative consequences, for exam- 
ple, in the scattering rates of neutrinos within the proto-neutron star matter. 
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Figure 1: The chemical potential /i for protons (solid lines) and neutrons 
(dashed lines) as a function of the proton abundance Y for T = 2 MeV and a 
set of fixed pressures 

P, < P < 0.19 MeV/ fm 3 . For each isobar both curves 
for the proton and the neutron chemical potentials are obtained, respectively. 
The Gibbs construction for a given isobar is shown in the figure inset (dotted 
line). 




Figure 2: Pressure P and free energy density J 7 in terms of the relative 
density n/n®, for T = 2 MeV and several isospin asymmetries w. Solid lines 
represent the equilibrium values for stable matter, the break in the slope of 
these curves indicates the end point of the transition. P and T evaluated for 
states out of equilibrium in the homogeneous one phase are shown by dashed 
curves. 




Figure 3: The binodal construction as explained in the text, for several 
temperatures. The curves in the (P, Y) plane enclose the zone where two 
independent phases coexist in thermodynamical equilibrium. 
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Figure 4: Zero sound dispersion V z — ui/qaXT — 1 MeV for several asymme- 
tries w of the nuclear medium. Solid lines correspond to modes propagating 
either in the binodal coexisting phase with the higher density for n < n , or 
in the stable one phase at supra normal densities. Conversely, dashed lines 
represent collective propagation within the unstable homogeneous phase, as 
explained in the text. The label isos (isov) indicates the iso-scalar (iso- vector) 
character, for damped (d) and undamped (u) stable modes. Curves which 
change its character from iso-vector to iso-scalar are labeled mix. For the 
unstable iso-scalar curve the ratio \fj\/q is plotted. The average Fermi veloc- 
ity V F is shown for the stable high density phase (dot line), and the unstable 
homogeneous phase (dash- dot- dot line). 
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